Denaturation of Circular DNA: Supercoils and Overtwist 
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The denaturation transition of circular DNA is studied within a Poland-Scheraga type approach, 
generalized to account for the fact that the total linking number (LK), which measures the number 
of windings of one strand around the other, is conserved. In the model the LK conservation is 
maintained by invoking both overtwisting and writhing (supercoiling) mechanisms. This generalizes 
previous studies which considered each mechanism separately. The phase diagram of the model is 
analyzed as a function of the temperature and the elastic constant k associated with the overtwisting 
energy for any given loop entropy exponent, c. As is the case where the two mechanisms apply 
separately, the model exhibits no denaturation transition for c < 2. For c > 2 and k = we find 
that the model exhibits a first order transition. The transition becomes of higher order for any 
K > 0. We also calculate the contribution of the two mechanisms separately in maintaining the 
conservation of the linking number and find that it is weakly dependent on the loop exponent c. 

PACS numbers; 87.15.Zg, 36.20.Ey 



I. INTRODUCTION 

The thermal denaturation of DNA, whereby the two 
strands of the molecule separate upon heating, has been 
thoroughly investigated both experimentally and theo- 
retically in the last half century. This process is relevant 
for experiments such as polymerase chain reaction (PGR) 
[l|, , and for biological processes such as those taking 
place within a thermophilic bacteria 0, Q • The fraction 
of bound base pairs vs. temperature (the melting curve) 
is measured by means of fluorescence and UV absorbtion 
methods. A typical melting curve of chains of the order 
of thousands of base pairs is composed of a sequence of 
discrete steps, interpreted as indicating a series of sharp, 
first-order phase transitions corresponding to the local 
melting of regions with successively increasing GC con- 
tent. 

A prototypical theoretical model for studying this 
phase transition is the Poland-Scheraga (PS) model 5]. 
In this model, and for the case of a homopolymer DNA, 
the molecule is represented by an alternating sequence 
of rigid bound segments and flexible denatured loops. 
Their contribution to the partition function are energetic 
and entropic, respectively. The entropy S(l) of a loop of 
length I is of the form 



where A and s are constants and c is the loop expo- 
nent depending only on dimensionality and constraints 
imposed on the DNA chain such as excluded volume in- 
teractions. In the framework of the PS model, the nature 
of the transition is set by the value of c: for c < 1 no tran- 
sition takes place and melting is just a gradual process in 
which the fraction of bound base pairs is nonzero at all 
temperatures; for 1 < c < 2 the transition is of second 



order; for c > 2 it is of first order, i.e., the melting curve 
is discontinuous at the melting temperature Tc. It was 
shown relatively recently jGj that the excluded volume 
corrections in three dimensions yield c ~ 2.12. Therefore, 
the PS model predicts a first-order melting transition. 

The DNA is a double helix and in order to open a dena- 
tured loop the region in which it is embedded must be un- 
wound. This has no consequence for a linear DNA chain 
in thermal equilibrium, where the ends of the molecule 
are free to rotate. On the other hand, in circular DNA 
(such as plasmids) and in DNA with rotationally fixed 
ends the total linking number (LK), which measures the 
number of windings of one strand about the other, is con- 
served. In such cases the unwinding of one region must 
be compensated by the over-winding of another region. 

Two mechanisms have been suggested to absorb the 
extra linking number in the overwound regions: (a) in- 
creasing of twist ("Tw"), or overtwisting, in which the 
change in the average stacking angle accounts for the ex- 
tra windings 0,i, and (b) increasing of writhe ("Wr"), 
or supercoiling, where the backbone assumes a nonpla- 
nar shape that accommodates a nonzero LK [i-[l3. The 
Galugareanu- White-Fuller theorem implies that, during 
the melting process one has LK — Tw + Wr |14| - [l7| It 
has been shown that the two mechanisms have similar 
effects on the melting behavior: For c < 2, the melt- 
ing process becomes a smooth crossover with no phase 
transition. For c > 2, there is a phase transition of 
high order, where the singular part of the free energy 
scales with the reduced temperature t = {T — Tc) as 
Fstng ~ |t|('="^)/('="^). Thus the order of the transition 
diverges as c approaches 2 from above, and it becomes 
second order for c > 3. Unlike the full denaturation of 
DNA with free ends, the high-temperature phase here is 
composed of a critical fluid of microscopic loops coexist- 
ing with a single macroscopic loop [l3j . 
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FIG. 1: A typical configuration of tlie model 



In this paper, we investigate the general scenario where 
both mechanisms act simultaneously. We find that the 
nature of the transition is the same as that found in ei- 
ther of the two mechanisms separately. This observation 
is expected, although it is not guaranteed by the fact that 
the two limits (supercoiling or overtwisting alone) yield 
similar phase transition scenarios. To probe the inter- 
play between supercoils and overtwist we also calculate 
the linking number absorbed by overtwist at the critical 
point. 

The paper is organized as follows: In section 2 the 
model is defined and analyzed with a formalism some- 
what different from earlier accounts. In section 3 the 
results are presented, first for simplified cases and then 
for the full model. We then conclude in section 4 with a 
brief discussion of our results. 



no internal entropy although they do contribute to 
the overall entropy through their positional degree 
of freedom. 

• A loop of size I has an entropic contribution given 
by the Boltzmann weight ^{l) — Af^. Here s is a 
geometry dependent constant and A is a constant, 
usually termed the cooperativity parameter, reflect- 
ing both the normalization of the entropic contri- 
bution and the enthalpic cost of initiating a new 
loop. 

In addition, we assume an overtwisting elastic energy cost 
with an elastic constant n. 

The excess linking number residing on the double- 
stranded DNA segments (bound segments and super- 
coils) is calculated as follows: A unit length of a loop 
region increases LK by 1 , while a supercoil segment of the 
same length decrease LK by 1. The 1 : 1 ratio is assumed 
for the sake of simplicity, while considerations of univer- 
sality suggest that the results should remain qualitatively 
unaltered under different choices. Denoting by Lb, Lg and 
Li the total length of the bound segments, supercoils and 
loops, respectively, the excess LK in the double stranded 
regions is simply Li — Lg. This is compensated by an 
increase in the average stacking angle per unit length 
by A9 in the bound segments and supercoils combined, 
hence A0 = 0. Then, the elastic energy cost due 



Lb + L 

to overtwisting is k (Lb 
ing the Hamiltonian 



Lg){A9y 



yield- 



H — Lb^b + LgEg 



Lb -\- Ls 



Thus, for example, the Boltzmann weight of the config 
uration depicted in Fig IT] is given by 



II. THE MODEL 

In order to incorporate supercoils and overtwist, we ex- 
tend the PS model and assume that each configuration is 
composed of an alternating sequence of bound segments, 
loops, and supercoils, the latter being double stranded 
chains which carry writhe. It is assumed that supercoils 
form only within bound segments. A typical configura- 
tion is sketched in FigHJ The contribution of the three 
types of segments to the free energy can be computed 
using the following rules: 

• A bound segment of length I contributes to the in- 
ternal energy Eb < per unit length, and none to 
the entropy (due to the large persistence length of 
the double stranded DNA). Hence, the associated 



Boltzmann weight is e' 



with 13 = l/ksT. 



• A supercoil of length I contributes to the internal 
energy < per unit length with E^ > Eb, yield- 
ing a Boltzmann weight e^^^" = . It is assumed 
that like the bound segments, these segments carry 



7b ,s ,0 ,s ,b 



X e 



-/3k. 



(^l+'2~^l,l~'l,2~'2,l) 



Jb JS Jt) 



It is worth noting that some of the previously studied 
PS-type models can be formulated as special cases of the 
present model: Ls = corresponds to the case with over- 
twisting only [7 1; Lg = Li corresponds to the case with 
supercoils only ll|, [l^ and k = corresponds to a DNA 
with supercoils and no LK constraint 

The canonical partition function can now be written 

as 



Lb+L, + L, 



Zf^^o{Lb, Li,Ls)e 



(1) 



where Zi^^o{Lb, Li, Ls) is the partition sum with given 
Lb,i^s and k = (which is an ensemble more restricted 
than even the microcanonical ensemble, as there may be 
different Lb^i.s triplets that have the same energy). The 
correspondence with other models mentioned above can 
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be obtained from Eq.([T]) by taking the appropriate limits: 
The PS model is recovered when k = and v = (or 
= oo so that Ls = 0); taking v = 0, k > yields 
the partition sum for a model with overtwisting only ; 
finally, substitu ting n = oo, u > recovers the case with 
supercoils only 



A. Free energy 

We begin by calculating ^^=0(^6, Li.Lg). This is con- 
veniently done by first evaluating the grand canonical 
partition sum by means of a z-transform of 2'k=o. The 
canonical partition sum (expressed in terms of Li and 
Ls, or in terms of the fractions rrii = Li/L, i = b, s,l ) 
is then calculated using the inverse transform. Introduc- 
ing three fugacities, Zh,zi, and Zs, corresponding to the 
three length constraints, the resulting grand canonical 
partition function QK.=oizb, Zs, zi) can be expressed in a 
closed form as 

QK=o{zb,Zs,zi) = Zi^=o(Lb,Ls,Li) z{^''Zs'zf^' 

= l + Vizb,Zs)U{zi) + 

+V{zb,Zs)U{zi)V{zb,Zs)U{zi) 



1 



with 



U{zi) 

V izb,Zs) 

V{zb) 

W{Zs) 



l-Vizb,Zs)Uizi)' 



00 I I 



1=1 



V{zb) + V{zt.)W{zs)Vizb) 
l-Wizs)V{zb)' 



1=1 
00 

E 

1=1 



I I _ UJZb 



,IJ 



1 — LOZb 



1 - VZs 



(2) 

(3) 

(4) 
(5) 
(6) 



The functions U (z;) and V (zb, Zg) are the grand canon- 
ical sums of single-stranded (loops) and double-stranded 
(bound and supercoiled) segments, respectively. Simi- 
larly, V {zb) and W (zs) denote the grand sums for bound 
and supercoiled segments, separately. Here $c (9) is the 
polylogarithm function of order c, which is analytic ev- 
erywhere except for a branch-cut for q > 1. The behavior 
of this function at q — I depends on c: If c < 1, ^c{q) 
diverges as 9 — > 1^. If c > 1, ^c{q 1^) = Cc where 
is the Riemann zeta function [31 . The behavior of $c('?) 
near q — I determines the nature of the phase transition 
investigated here, as will be shown below. In deriving 
Eq.(I2]) we take Z'„=o(0,0,0) = 1 and assume that the 



chain contains at least one loop and one bounded seg- 
ment. This assumption simplifies the numerator of the 
resulting expression in ([2]) and it has no effect on the 
resulting thermodynamic properties of the model. 

The canonical partition function is found by inverting 
the z-transform using a Cauchy integral: 



Zf^^a{Lb,Li,Ls) — - — : 

' ZTTl 



Qi^=o{zb,Zs,zi) 

Zj, Zs 



dzbdzsdz^ 
(7) 

All integration contours encircle the origin and contain 
no other singularities. Using Eqs.([2]|6|) we find 



^o{zb,Zs,zi) 



1 

UJZb 



1 



1 — vz, 

Qk=o has a simple pole in Zb set by 
1 1 

UJZb 1 — l^Zs 

yielding 



A$c {sZl) 



A^c {szi) = 0, 



1 



1 



1 - VZs 



-4$, [szi] 



(8) 



(9) 



Note that Eq.® is equivalent to Eq.(2) in [Tl|. The z& 
contour in Eq.([7]) can be deformed to encircle the pole 
given by Eq.®, yielding 



Zk=o {mb,mi) = 



with 



2m 



(10) 



Ff,^o{zs,zi,mb,ms) = m;, log [z^ (z^, z;)] + log (z^) 
-f (1 -mb-ms) log (z/) (11) 

up to logarithmic corrections in L. Here we used the fact 
that nib + ma + mi = 1. In the thermodynamic limit the 
integral in Eq. (flO|) can be evaluated by considering the 
saddle point of i^K=o with respect to z/ and Zs, 



dP^^o _ mbu/{l - vzs) TTh 

dzs 1 + (1 - l'Zs)A<^c{sZl) Zs 

OFk^o _ mb(l - iyZs)A^c^i{szi) 
dzi ~ zi[l + {l- i^Zs)A<^>^{szi)] 

1 — mb — ms 
zi 



(12) 



(13) 



where we used the identity -^'^c (<z) = \'^c-i (q)- After 
some algebra Eq. (fT2| yields 



l + x - y (l + x) + Axri 



where we define 



mb 



, r) = A<^ciszi). 



(14) 



(15) 
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It can be seen (as 1,77 > 0) that z* < v~^, therefore the 
Zs integration contour can be deformed to pass through 
this saddle point without encirchng the singularity of F 
at = i/-^ Eq.dni) yields 



1 



A^c-i{szi) 
\-^+A^,{szi) 



1 



1 



(16) 



The LHS of Eq. ([T6|) is monotonically increasing with zi 
(see Appendix A). For c < 2 this equation has a solu- 
tion, Zi, for any value of mjj and x due to the fact that 
<^c-i{szi) diverges at szi — 1. However, for c > 2, the 
LHS reaches a finite value for zi — s^^ (the branch point 
of ^c{szi)). Therefore, for a given value of x, there is no 
saddle point for values of Wfc below a critical threshold 



given by 



ml \x) = 



1 + 



(17) 



Here z* is obtained using Eq. p^ with 77 replaced by AC,^- 
For values < m!f' (x) the zi integral is equal to the 
value of the integrand at the branch point zi — s~^, as in 
the canonical treatment of the PS model [19]. The inte- 
gration procedure involves more than simply evaluating 
the integrand at the singularity closest to the origin. De- 
tails are given in Appendix B. 

To calculate the canonical partition function Z{L) 
given in Eq.([T]) the overtwist term should be added to 
the free energy, yielding 

F {zs,zi,mb,ms) = mb log [z^ (z^, z;)] 

-I- ms log (zs) + (l-mb- rus) log (z;) 



13k 



(1 - Wb - 2ms) 



(18) 



rub + rus 

This full free energy needs to be minimized with respect 
to all of its arguments. The minimization with respect 
to Zs and zi is the same as for -F^^o a-nd the results are 
given by Ea. (ll4ll6p . The minimization with respect to 
TOb and rus is discussed in the next section. 

In summary, the fugacities Zg , z; in the thermodynamic 
limit (denoted by z* and Z;*) as functions of the bound 
and supcrcoiled segment fractions, mb and mj, are given 
by Eqs. (|14ll6p . Hence we can express the Landau free 
energy, Eq. (fT5)l . as a function of the densities rub and 
TOs only. In what follows it will be occasionally more 
convenient to express the dependence on TOs through the 
fraction x = mb/ms as defined in Eq. (|15p . Then the 
Landau free energy can be written as 



F{mb,ms) = Fi^^Q{mb,x) + 13 k 



(l-mb - 2ms) 



mb -I- ms 



(19) 



Fn^a{mb,x) = Ff,^o{mb,ms,zi{mb,x) ,Zs{mb,x)) 
= rub log [zb{zs, Z()] H log (zs (x, z/)) 

X 

+ (1 -mb^) log (z,), (20) 



where Zb is given by Eq.®, Zg is given by Eq. p4l) and 
z/ = zi{mb,x) is given by Eq. (fTG)) . Note that unlike 
other Landau free energies which are analytic in the order 
parameter, here F„=o(wb,a;) is a non-analytic function 

of rUb and x along the line defined by m['^^ (x) . We shall 
denote this line of singularities by F. In the mb — mg 
plane the expression for F is 



F : {mb,ms) 



mi^hx) 



ix) 



x e (0,00). (21) 



For points to the left of F (i.e. those for which nib < 
mlf^ (x) so that z; = s^^) and for a given x the free 
energy is linear in rub, while for points above F it has a 
more complicated form, hence Ff^^olmb^ms) is singular 
along F. Below we will explore this non-analyticity in 
more detail and show that it is closely related to the 
non-analytic behavior of the free energy as a function of 
temperature at the transition point. 



III. RESULTS 

After introducing the Landau free energy F(mb,ms) 
and arguing that it is singular on a line (F) in the 
{mb,ms) plane, we move on to study the nature of the 
phase transition for different values of k. Three cases are 
of interest: 

• K ~ Q: Here overtwisting has no cost and the chain 
is equivalent to a linear chain with supercoils freely 
spread within, with no linking number constraint. 
We will see that in this case the transition (which 
exists only for c > 2) is first order as in the standard 
PS model. 

• K 00: Here overtwisting is forbidden. This is 
the case with supercoils only which was analyzed 
in Ref. [l3[ and found to exhibit a continuous tran- 
sition of order with a singularity in the free 

energy which scales as ^ t'^'^^^) / ('^~'^) _ Here we will 
outline the derivation of this result within the cur- 
rent approach. 

• < K < cxi: In this case supercoiling and overtwist- 
ing coexist, yielding a different free energy mini- 
mum. Yet, it is shown that the nature of the tran- 
sition remains the same as for k — 00. 

To quantify the interplay between supercoiling and over- 
twisting, we calculate the fraction of the linking number 
accommodated by overtwist at the transition point 



r{K) = mic - msc = 1 - mbc 



2m, 



(22) 



Clearly, r{K — co) = as no overtwist is allowed in this 
limit. Below we derive an explicit formula for r(K — 0) 
and calculate r(0 < /t < 00) numerically. 

Throughout the paper the parameters used in the fig- 
ures are Eb = -3, E^ = -2, s = 5, A = 0.1. 
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A. K = 

The densities nis and mf, are found by minimizing 
F {mi,, x) = i^K=o {'Ttbi x)- This is equivalent to minimiz- 
ing F (zg, zi, mi,, TTig) with respect to all of its arguments, 
as F is obtained from F by minimizing it with respect to 
the fugacities Zs and z;. Using Eq- ipO)) and minimizing 
F (mf), ms) yields 



Zs {mi,, x) = zi {mi,, x) = Zf, {mi,, x) . 



(23) 



i.e., the system is described by a single fugacity. This is 
expected, since for k = the grand canonical partition 
function of the full model could have been derived with 
a single fugacity corresponding to the single constraint 
Lb + Ls + Li = L. 

Substituting Eq.dSS]) in Eq.dSO]) one finds 

^^(r)-log [zi (m,(r),x(T))], 

where mi,, ms and zi can be calculated by solving 
()9ll4ll6p using Eq. ([23|) . Hence the non-analytic behavior 
of zi results in a singularity in F. As mentioned above, 
for c < 2, zi {mi,,ms) is an analytic function, therefore 
there is no phase transition in the system. For c > 2, zi 
increases monotonically from zero with temperature as 
long asT <Tc, where using Eq.® the critical tempera- 
ture Tc is given by 



1 



n -1 



1 - V (Te) /S 



For T > Tc, zi = is a constant and hence F{mi,, x) is 
constant, independent of T. As T — Tc from below, the 
free energy approaches the transition point with a non- 
zero slope, i.e., the transition is first order. This can be 
seen by differentiating Eq.® with respect to T ~ log(iu;) 
using relation (1^^ . 

The equilibrium values of and ms as a function of 
T define a trajectory in the mi, — ms plane, as depicted 
in Figl2] The starting point of this trajectory at T = 
is to the right of the singular line F since mt > m^^^ {x) . 
As T increases, mi, decreases. At T = Tc where the 
trajectory intersects F, the singular line defined above, a 
phase transition takes place. 

The fact that the transition is first order can be veri- 
fied as follows: the intersection of the trajectory with F 
takes place at a certain Xc — mhcjmsc where m^c and 
TOsc are the critical fractions of the bound segments and 
the supercoils on F, respectively. Since the minimum 
of the Landau free energy at the critical temperature is 
obtained at {mi,c,msc), the slope of this free energy van- 
ishes in all directions, i.e. dF {mi,c,Xc) — 0. As stated 
above and can been seen by inspecting Eq. ipn]) . for points 
to the left of F, where m;, < mi,c and T > Tc, the free 
energy F {mi,,x) is linear in mi, for fixed x. Hence the 
slope of F{mi,,Xc) for mi, < mi,c must be 0. This im- 
plies, in particular, that F {mtc, Xc) = F{0,Xc), namely 
a phase coexistence between bound and unbound phases. 




0.4 0.6 0.8 
m. 



FIG. 2: (Color online) The {mi,,ms) trajectories as a function 
of temperature of for k — (green full and dotted line), k — 1 
(red dashed line) and n — oo (cyan dash-dotted line). For all 
lines c = 2.5. At T = the chain is totally bound, so mt — 1 
and rus — for all trajectories. As T increases nit decreases, 
intersecting at T = Tc the singular line F (Eq. (|21|) '). drawn 
above as a thick black line. On the n — trajectory there is a 
coexistence between a bound phase and a denaturated phase 
at T = Tc, along the dotted green line, while for T > T the 
system is effectively unbound with mi, = = 0. For k = 1 
(and any k > 0) the trajectory continues smoothly across the 
the singular line and reaches nit — rUs — only at T = oo. 
For n — oo the trajectory is linear due to the simple relation 



nis = i(l 



Note that there is no free energy barrier between the two 
phases. As depicted in FiglSlJa) and can be verified by 
Eq. (f20|) above Tc the slope of F{mi,,Xc) for mf, < mi,c 
is positive and hence mf, = ms = is the minimal so- 
lution, so the system is in the unbound phase. As will 
be discussed below, setting k > eliminates the coex- 
istence and yields a unique free energy minimum at all 
temperatures (see Fig|3Ub)). 

Let us now consider the overtwist linking number 
r{K — 0). The value of a; = miy/ms at criticality can 
be calculated using Eqs. (|12l23|) : 



-7(1 



l + {l-vs-^)ACc' 



Solving Eq. pTj) and msc — m^cXc ^ for mf,c and TOs 
obtain the overtwist linking number at Tc as 



we 



r{K = 0) = 



ACc 



1 



1 + 



-1)' {ACc + AQ-i) 



Depending on parameters in this expression, r(K — 0) 
can be either positive or negative. Specifically, for the 
parameters used in Figl2]the value is r(K) — —0.26343 < 
0, implying that, the length of the supercoiled regions at 
the phase transition point exceeds the length needed to 
compensate for the linking number released by the loops, 
resulting in undertwisted bound segments. 
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FIG. 3: The Landau free energy F{mb,ms) as defined in 
Ea. p9|l along the line nib ~ Xcms where Xc = for c = 2.5 
and (a) ft = ; (b) k = 0.01. While the free energy at the 
phase transition point for k = has a continuum of min- 
ima in the interval mt < mj'^'(a;), for k — 0.01 (and for any 
K > 0) there is a unique minimum of the free energy at all 
temperatures. 



FIG. 4: nit vs. T for various c and k. The signature of 
the second order transition is the non differentiability of the 
curves at T = Tc when c > 3. For 2 < c < 3 the melting 
curve is smooth with a higher order singularity at T = Tc. It 
can also be seen that the transition sharpens as k decreases. 



decreases as c increases and yields a 2"'' order transition 
for c > 3. This can be seen by expanding Eas. (|12ll3l24p 
near the critical temperature, where sz; = 1. Setting t = 
Tc-T, Srub = mb-mbc, Szi = zj-s"^ and Szs = Zs-z^c, 
where mbc and the values of and Zg at T^, and 

using the identity $c(l — ^) ~ Cc ~ S'^^^ yields below the 
critical temperature [t > 0) 



6mb '■'^ Szs ^ Szf 



Hence Szi ~ t"-^. Expanding Eq. (|25|) to appropriate 
order in t and Szi yields 



B. 



When K — oo, overtwisting is forbidden. The con- 
servation of the linking number now implies nis — mi, 
therefore = ^— Minimizing Ea. ([T5)) with respect 
to rub and using the linking number constraint yields the 
relation 



F = F{Tc) + at + /3Sz^ 



o e 



Zb = 



(24) 



which implies that two fugacities are needed, accounting 
for the two constraints on the linking number and the 
total chain length. Indeed, in previous accounts of this 
model the derivation was conducted using two fugacities 
[ni[il- Inserting Eq.dH]) into Eq.((ni) yields 

F = i log (z,) -f i log {zi) = log {zb) . (25) 

In Ref . [l^ this case was analyzed and the transition was 



found to be of order 



which diverges as c 2+, 



F {Tc) + at + pt— + O if 



(26) 



Above the critical temperature [t < 0) 6zi = and hence 
the -th derivative of F diverges as temperature ap- 

proaches Tc, constituting a phase transition of the same 
order. The fact that the transition becomes more pro- 
nounced (of lower order) as the loop exponent c increases 
can be appreciated by inspecting Fig|3] which shows nib 
as function of temperature for c below and above 3. 



C. < K < oo 

In this case, as for k = 0, both mf, and mg are set 
by minimizing the Landau free energy given in Eg. (1191) . 
Here, however, there is no simple relation between the 



7 



fugacities 



= 



= 



dF 
drub 

log 

dF 
drUg 

-13k 



13k 



(1 - ma - 2ms) (1 + rrib) 
{rrib + rusf 



= log 



(1 — TOb — 2ms) (1 + 3mf, + 2tos) 



(27) 



(28) 



(mf, + rris) 

These equations, together with Eqs. (|14|16p for T < Tc, 
and Eqs.(fT4|) and zi — for T > Tc , set the value 
of the order parameter in the thermodynamic limit. 
Inserting Eas. (|27l28p into Eq-dUl) yields 



(1 - nib - 2tos) 
(mfc + nisf 

- oo shows that here, 
-1 



i^(T) = log [zfc (mfc,m,)]-/3K 

Repeating argument used for k - 
too, the order of the transition is 

We observe in Figl2]that the trajectories for k = 0, 1, cx) 
in the (mb, vis) plane intersect at a single point (mj, m*). 
In fact, this special point is common to all such trajec- 
tories with arbitrary k: Let T* be the temperature for 
which the minimum of the free energy F (mbjnis) sat- 
isfies m* = ™'' for some k. Then, for any other k, 
the minimum of the free energy at T* is also given by 
(771^,771*), because the K-dependent part of the free en- 
(i-mb-2ms) yg^jjjgj^es (and hence is minimal) when 



ergy k 



We now consider the overtwist linking number at crit- 
icality for < k < oo. This number, r = 1 — rubc — 'inisc, 
cannot be obtained analytically. In FiglS]we present the 
numerically calculated r(K)/r(0) ratio for two values of 
c. r{K) depends weakly on c and could be either positive 
or negative, depending on the parameters of the model. 
However, for a given set of parameters, the sign of r{K) 
does not change with k. In order to demonstrate this 
point, consider the special point in Fig [2] where all tra- 
jectories for different k intersect at a shared temperature 
T* and note that this point is the borderline between neg- 
ative and positive r{K) on each trajectory. Therefore, if 
the parameters are such that the intersection is to the left 
of the singular line F, then for a given k the critical tem- 
perature satisfies Tc < T* and hence r(K) < 0. If, on the 
other hand, the intersection is to the right of the singular 
line then r{K) > for the same reason. In addition, if the 
parameters are such that the intersection is to the left of 
the singular line, a corollary follows that T* = yj'"- o)^ 
which in turn implies Tc'^^^'^ < Tc'^~^\ Recalling that 
K = refers to the case with no LK conservation, this 
is in agreement with the experimental evidence that im- 
posing circular topology reduces the melting temperature 
There is no equivalent statement in the other case 
in which the intersection is to the right of the singular 
line. 



0.8 
o 0.6 

'5" 

^0.4 
0.2 



c=2.5 

---c=10 



12 3 4 
K 

FIG. 5: r(«:)/r(0), where r = 1 — nitc 
values of c. It can be seen that r(K)/r(0) depends only weakly 
on c, and that, as expected, it decays monotonically from 1 
for K = to at K = oo. 



IV. CONCLUSIONS 

In this paper we analyzed the thermal denaturation of 
a circular DNA molecule, in which the linking number is 
conserved. Within the framework of the Poland-Scheraga 
model, we have considered the two possible mechanisms 
for conserving the LK: writhing (forming supercoils) and 
(over)twisting.The denaturation transition is studied for 
arbitrary values of the elastic constant k associated with 
the overtwist elastic energy. We found that the model 
exhibits no transition for c < 2 and a high-order, contin- 
uous transition for c > 2, k > 0. The singular part of the 
free energy was found to scale as t'=-'^ with t — Tc — T, 
yielding a transition of order . The order of the 

transition diverges as c approaches 2 from above, it de- 
creases with increasing c and it becomes second order for 
c > 3. The model with k = behaves differently, exhibit- 
ing no transition for c < 2 and a first order transition for 
c > 2. Similar observations were reported before for the 
limiting cases restricted to supercoils only Jl3] and over- 
twist only [20| . 

The canonical analysis carried out here brings new in- 
sights. For example, the first-order transition which takes 
place for k = and c > 2 is found to be rather special 
in that it does not have a metastable region (see Fig|3K). 
This is true also for the original PS model. In addition, 
the analysis of the {mb,ms) trajectories unveiled a k- 



independent special point (777^ , 

n(K>0) 



1*) which in return led 
to the prediction that Tc'"^"' < j'j''"'^^ (the melting tem- 
perature reduced by circular topology) for a wide range 
of parameters, in line with an earlier experimental obser- 
vation. 

The model considered in this paper corresponds to a 
homogeneous circular DNA chain, while biological DNA 
molecules are heterogeneous. However, through the Har- 
ris criterion (2l| we find that the disorder is irrelevant 
for K > and c < 3, where the specific heat exponent 
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a = 2 — i = ^— I is negative. Therefore we do not 
expect the sequence heterogeneity to change the nature 
of the phase transition and the associated critical ex- 
ponents. As the actual value of the loop exponent was 
estimated to be c « 2.12 Q, our analysis should be valid 
for sufficiently long, real DNA chains. 

Previous accounts on denaturation of circular DNA 
have found that a macroscopic loop is formed above Tc, 
reminiscent of Bose- Einstein condensation. Although not 
discussed here, we expect a similar phenomenon in the 
combined model of supercoils and overtwist, and it would 
be be interesting to analyze the linking number exchange 
between the macroscopic loop, the microscopic loops, and 
the supercoiled and overtwisted segments. 
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V. APPENDICES 
A. Appendix A: Monotonicity of the RHS of 

Eq.dieD 



We wish to show that 



A$c-l(sZ/) 



is a monotonically increasing function for zi e (0, s ^). 
Using Eq.(fT4|) we can write 

5 ('7,2;) 



\J{l + xf + Ax-q - (1-1- x) 



where r] = A$c {szi) increases and g (r], x) decreases with 
zi. Hence we can write 



1 <^c-l{sZl) 

1-1-5(77,2;) '^c{szi) 



where the first factor is an increasing function of z/. It 
is thus sufficient to show that the second factor also in- 
creases with zi. To this end, we differentiate this term: 



d 


'^c-i{szi)' 




''i>c-2{sZl)^c{sZl) 


dzi 


'^c{sZl) 




*c(sZ/)2 



^c-l{sZlf 



^c{sZlY 

Now we show that the numerator of the derivative, de- 
noted by S {zi) is positive, by expressing the poly loga- 
rithm function explicitly as a power series of the variable 



y = szf. 

00 

S(^0 = E 



yl yk yl yk 



k,l=l 
00 / h. 

— y y 



lc-2 l^c lc-1 

00 



k<l 
°° „/+fe 



11 11 



k<l 



This demonstrates that f{zi) is a monotonically increas- 
ing function of z;. 



B. Appendix B: Branch-cut integration 

We wish to evaluate the integral for the partition func- 
tion with K = given in Eq. pU]) 

2m J 

with Zs {x, zi) given by Ea. (ll4p . Li satisfy L11 + Li + Ls = 
L, Lb/ Ls — X and rrii = Li/ L (i = b, s). Defining y = szi 
yields 



ip {Lb, x) 



s ' 
27ri 

2m ' 



A<^c [y] 



-dy 



Eq. ([T7|) defines the value of nib = Lb/L below which 
the integrand has no saddle point, therefore the integral 
should be evaluated by another method. The integration 
contour can be deformed to the contour depicted in Fig|6l 
composed of the following segments: 

(/) : [R-ie,l-ie] 



(///) : [l + ie,R + ie] 
{IV) : {i?e"' :6<e<2Tr-S} 



(29) 



where R ~> 00 and tg {6) ~ . We wish to show now 
that the only contribution comes from the vicinity of the 
branch point: To sec that the contribution of {IV) is 
negligible we note that for \y\ — > 00, \^c{y)\ r(c+i) 
[2^ . From Eq. p^ we see that for |y| — >■ 00 z^ -H- 1 and 
^ ^/^c{y) so that for large enough L and a; < 00, 

lim r/(i?e^^)]^~ hm ^""^^^^H" < lim i- . 
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CI 



(III) 




(I) 



FIG. 6: The original contour of integration CI can be de- 
formed to the contour which is given by Eq. (|29|) and is de- 
scribed in the text of Appendix B. The branchpoint y = 1 is 
marked here by a filled circle. 



Hence segment {IV) of the contour has no contribution. 
Along segment (//) the function I{y) is analytic and 
hence the integral is of order e and can be taken to be 
arbitrarily small. 

I{y) can be written as a power series with only real 
coefficients, and hence I (y*) — I* (y) , where y* is 



the complex conjugate of z. 
(/) + {III) yields 



fntegrating along segments 



I- / 



Im 



i{yr 



dy. 



Defining I{y) = A{y)e''^^y\ where A{y) and ^(y) are real 
functions, we can write /m [/(y)^] =A{y)^sin[L'${y)]. 
As ^'(y) is a smooth function for y > 1 in the thermody- 
namic limit the oscillations in the sin() function average 



out to zero. Therefore, the only contribution to the inte- 
gral comes from the vicinity of y — 1, where ^'(y) has a 
discontinuity in some derivative. 

The function I{y) has a pole at y = where I{y — >■ 
0+) — >■ 00, so that I'{y) < near the origin. When 
c > 2, the fact that there is no saddle point for < y < 1 
implies that /'(I) < as well, and hence A'(l) < 0. 
For Sy = y — I ^ I the imaginary part of the poly- 
logarithm function is approximately Im [$c(l + Sy)] « 
aSy'^"^ with a = n/r{c), and therefore ^'(1 -I- Sy) 
aSy'^^^ with a — (a/A(l)) x dl/d^c{y)- Combining 
these observations yields to leading order /(I + Sy) « 
A(l) exp [-bSy + ihSy"-^] with b = -A'(l)/A(l) > 0. 



« A(l)^ e-^^^y sm{LhSy'-^)d5y. 

The contribution to the integral comes from a region of 
size Sy ^ j^, so that the upper limit can be stretched to 
oo without affecting the result. Rescaling by y = hLSy 
we obtain 



^ « '-!lK{l)HhL)-' j\-ysn.y^^]dy. 



As c > 2, in the thermodynamic limit L — > oo the ar- 
gument of the sin() function is small and it can be ex- 
panded. 



^ « s^'A(l)^ / e-yy'-^dy 

Jo 

d {bL)^^'^ r (c) r . , , r 

= — ^ — '- ^s^'A 1)^ . 

ttL 

Substituting back zi — y/s yields 
ttL 

where z^ = Zg {zi,x) is given in Ea. ((H|) . Hence up to log- 
arithmic corrections the free energy is given by its value 
at the branch-point. 
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